Nature of the ordering of the three-dimensional XY spin glass 
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Spin and chirality orderings of a three-dimensional XY spin glass are studied by extensive Monte 
Carlo simulations. By calculating an appropriately defined spin-overlap distribution function, wc 
show that the finite-temperature chiral-glass transition does not accompany the standard spin-glass 
order, giving support to the spin-chirality decoupling pictureD Critical behavior of the chiral-glass 
transition turns out to be different from that of the Ising spin glass. The chiral-glass ordered state 
exhibits a one-step-like peculiar replica-symmetry breaking. 
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Considerable attention has recently been paied to the 
ordering of the XY spin-glass (SG) model with two- 
component vector spins. The reason of such interest 
is probably twofold. First, the XY SG has found ex- 
perimental realizations, not only in SG magnets with 
an easy-plane-type anisotropy, but also in ceramic high- 
T c superconductors with the d-wave pairing symmetry 
which can be regarded as a random Jesephson network 
jij. Recently, the latter system has been studied exten- 
sively H-§|], and the XY SG is expected to serve as a 
reference model in interpreting the experimental data. 

The second reason of interest in this model is more con- 
ceptual. Namely, this model is the simplest realization of 
the random and frustrated models with vector internal 
(or spin) degrees of freedom. Ever since the pioneering 
work of Villain ||, this type of model has been known 
to sustain nontrivial "chiral" degrees of freedom corre- 
sponding to the sense of the noncollinear ordered-state 
structure stabilized by frustration. The ongoing contro- 
versy mainly concerns with the manner how the spin and 
the chirality order in such chiral systems. 

Earlier numerical studies suggested that the XY SG 
in less than four dimensions did not exhibit any finite- 
temperature transition J^,|| . In a series of papers on the 
XY SG in two and three dimensions, however, Kawa- 
mura and Tanemura observed a novel possibility argu- 
ing that the chirality associated with Zi spin-reflection 
was "decoupled" from the spin assocaited with 5*0(2) 
spin-rotation on sufficient long length and time scales 
(spin-chirality decoupling) ]^-pd|] . More specifically, they 
claimed that, in two dimensions (2D), while both the spin 
and the chirality order simultaneously at zero tempera- 
ture, the associated spin and chirality correlation-length 
exponents are mutually different, i.e., v s ~ 1 for the spin 
and v K ~ 2 for the chirality |9pC[|. In 3D, they suggested 
the occurrence of a novel chiral-glass transition at a fi- 
nite temperature, where only the chirality exhibited a 
glassy long-range order (LRO) without the conventional 
SG order g§g. 

For the 2D XY SG, general concept of such a spin- 
chirality decoupling was recently challenged. Koster- 
litz and Akino claimed on the basis of their numerical 
domain-wall renormalization-group (DWRG) calculation 
that the spin- and chiral-correlation-length exponents at 



the T = transition are common, i.e., v s = v K ~ 2.7 
]l2[ , while Ney-Nifle and Hilhorst made an analytical ar- 
gument for a certain 2D XY model that the equality 
v s = v K should hold By contrast, direct Monte Carlo 
simulations on the 2D XY SG have invariably suggested 
v K ~ 2 > v s ~ 1, apparently supporting the spin-chirality 
decoupling picture |l4]-|l7|]. 

In 3D, Granato recently suggested on the basis of a 
dynamical simulation of the ± J XY model that the spin 
and the chirality order simultaneouly at T ~ 0.4J with 
v s = 1-2(4) |Q. This observation might indicate the ab- 
sence of the spin-chirality decoupling in 3D. Meanwhile, 
on the basis of a numerical DWRG calculation, Maucourt 
and Grempell suggested that the SG order might occur at 
a nonzero temperature below the chiral-glass transition 
temperature, i.e., < T$q < Tqg P^l - 

The purpose of the present Letter is to make further 
extensive Monte Carlo (MC) simulations on the 3D XY 
SG to clarify some of the issues concerning this model. 
Our main goal here is twofold. First, we wish to study 
the controvertial issue mentioned above, i.e., whether the 
chiral-glass order accompanies the standard SG order or 
not. Second, we study whether the ordered-state of the 
3D XY SG, either the chiral-glass or the spin-glass, ac- 
companies the replica-symmetry breaking (RSB), and if 
so, to reveal its nature. We have found a strong numerical 
evidence that the chiral-glass transition does not accom- 
pany the standard SG order, and that the chiral-glass 
state exhibits a one-step-like peculiar RSB. 

The model we consider is the 3D XY (plane rotator) 
model, defined by the Hamiltonian 



h — *y ' Si 

<ij> 



cos 



ty, (1) 



where Si is the two-component unit vector at the i-th site 
on a 3D simple cubic lattice, while the nearest-neighbor 
random coupling is assumed to take either the value J 
or — J with equal probability (±J distribution). The lo- 
cal chirality may be defined at each elementary plaquette 
a of the lattice by, 



n Q = (1/2V2) sgn(Jy)sm(0i - 9 3 ), 



(2) 
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where the summation is taken over four bonds surround- 
ing the plaquette a in a clockwise direction. The chirality 
is a pseudoscalar invariant under the SO (2) spin rotation, 
but changes its sign under Z 2 spin reflection. 

Using the temperature-exchange MC method [p0| , we 
have performed a large-scale MC simulation superceding 
the previous simulations, in that we have succeeded in 
equilibrating the system down to the temperature con- 
siderably lower than those attained in the previous ones 
fTT| . By running two independent sequences of systems 
replica 1 and 2) in parallel, we compute a scalar chiral 
overlap q K between the chiralities of the two replicas by 
1 K = Jn K a K a , as well as a spin-overlap tensor q^ v 
between the /i and v (/i, v=x,y) components of the spin 

by q^v = 7f J2i SipSj®. Then, in terms of these over- 
laps, we calculate the Binder ratios of the chirality g Kl 
and of the XY spin g s , defined in the standard man- 
ner: See Ref. |ll],|2lj for detailed definition. The lattice 
sizes studied arc L = 6,8,10,12 and 16 with periodic 
boundary conditions. Sample average is taken over 1500 
{L = 6), 1200 (L = 8), 640 (L = 10), 296 {L = 12) and 
136 (L = 16) independent bond realizations. 

As can be seen from Fig. 1(a), the Binder ratio of the 
chirality g K exhibits a negative dip which, with increas- 
ing L, tends to deepen and shift toward lower tempera- 
ture. Furthermore, g K of various L cross at a temperature 
slightly above the dip temperature Tdip on the negative 
side of g K , eventually merging at temperatures lower than 
Tdip- We note that the observed behavior of g K is similar 
to the one recently observed in the 3D Heisenberg SG 
plfl . As argued in Ref. pi] ], the persistence of a nega- 
tive dip and the crossing occurring at g K < is strongly 
suggestive of a finite-temperature chiral-glass transition 
at which g K (T^ G ) and <? k (Tcg) take negative values in 
the L — ► oo limit. In the inset of Fig. 1(a), we plot the 
negative-dip temperature Idi p (L) versus 1/L. The data 
lie on a straight line fairly well, and its extrapolation 
to 1/L = gives an estimate of the bulk chiral-glass 
transition temperature Tqq/J ~ 0.41. (More precisely, 
Tcg(L) should scale with L 1 ^" where v K is the chiral- 
glass correlation-length exponent. As shown below, our 
estimate of v K ~ 1.2 comes close to unity, more or less 
justifying the linear extrapolation employed here. Ex- 
trapolation with respect to L 1 / 1 - 2 yields TqgI J ~ 0.38.) 
Our present estimate of Tqg is somewhat higher than 
the previous estimate of Ref. ]ll{ ], but is in agreement 
with the estimate of Ref. jig] . In Ref. |ll| , Tqg was de- 
termined as a point where g K appeared to merge on the 
positive side of g K , yielding an estimate Tqg = 0.32(3). 
However, since the merging of g K (L) develops very slowly 
with L, it is not easy to precisely locate the merging point 
and we believe our present estimate of Tqg is more reli- 
able than that of Ref. Q . 

In sharp contrast to g K , the Binder ratio of the XY spin 
g s decreases monotonically toward zero with increasing 
L, without a negative dip nor a crossing, suggesting that 
XY spin remains disordered even below Tcg- 
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FIG. 1. The temperature and size dependence of the 
Binder ratios of the chirality (a), and of the spin (b). Inset of 
Fig. (a) displays the the nagative-dip temperature versus 1/L. 
The solid and broken lines are the best fits assuming the 1/L 
and 1/L 1 / 1 ' 2 dependence, respectively. 

We also calculate the chirality autocorrelation function 
defined by 



(3) 



where < • • • > and [• • •] represent the thermal average 
and the sample average, respectively. MC simulation 
is performed here according to the standard Metropo- 
lis updating. The starting spin configuration at t = to is 
taken from the equilibrium spin configurations generated 
in our exchange MC runs. The result, shown in Fig. 2 
on a log-log plot, suggests that the chiral-glass transition 
occurs at T/J — 0.39(3), in agreement with the above 
estimate. Below T — Tcg> C K {t) shows an upward cur- 
vature indicating that the chiral-glass state has a rigid 
LRO characterized by a finite chiral Edward-Anderson 
order parameter q^Q > 0. In order to see the possible 
finite-size effect in our estimate of Tcg , we also take lim- 
ited amount of data for L = 20 (30 samples only) for 
comparison, and have checked that the present estimate 
of Tcg is indeed stable. 
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FIG. 2. Log-log plot of the time dependence of the equilib- 
rium chirality autocorrelation function for L — 16 at several 
temperatures. The best straight-line fit of the data, repre- 
sented by the solid line, is obtained at T/J ~ 0.39 with a 
slope ~ 0.077. 

With setting Tqq/J = 0.39 as determined above, we 
perform the standard finite-size scaling of the chiral-glass 
order parameter [< q 2 K >] and of the chiral autocorrela- 
tion function, to estimate various chiral-glass exponents. 
We then get v K = 1.2(2), r) K = 0.15(20), z K = 7.4(10). 
As compared with the estimates of Ref. v K ~ 1.5 
and rj K ~ —0.4, n K is considerably larger, mostly due to 
the difference in the estimated Tqg values. We note that 
the present u K and r\ K values are clearly different from 
those of the standard 3D Ising SG, suggesting that the 
chiral-glass transition lies in a universality class different 
from the Ising class. 

The observed deviation from the Ising behaviors, not 
only in the critical properties but also in g K , is likely to 
arise from the long-range nature of the chirality-chirality 
interaction ||. Indeed, in case of the 3D XY SG coupled 
to fluctuating gauge fields where the chirality-chirality in- 
teraction becomes short-ranged due to screening, the ex- 
ponents turn out to be close to the standard 3D Ising 
SG values and the negative dip in g K is absent, in sharp 
contrast to the present results p2| . 

In order to probe the possible RSB in the chiral- 
glass ordered state, we display in Fig. 3(a) the distribu- 
tion function of the chiral-overlap defined by P{q' K ) = 
[{S(q K ~ q' K ))} at J/T = 2.95 well below T CG - The ex- 
istence of a growing "central peak" at q K = for larger 
L, in addition to the standard "side-peaks" correspond- 
ing to ±9cg, suggests the occurrence of a one-step- like 
peculiar RSB in the chiral-glass ordered state. Similar 
behavior was recently observed in the chiral-glass state 
of the 3D Heisenberg SG jilj . The existence of a nega- 
tive dip in the Binder ratio g K is fully consistent with the 
occurrence of such a one-step-like RSB Q . 

We next turn to the spin-overlap distribution. While 
the spin-overlap distribution is originally defined in the 
four-component tensor space, we introduce here the pro- 
jected distribution function defined in terms of the "diag- 



onal" overlap which is the trace (diagonal sum) of q^s, 

£w = ~E3 1, -3 9) . (4) 



<Zdiag 



The distribution function P(gdiag) is symmetric with re- 
spect to gdiag = 0. In the high-temperature phase, each 
q^y is expected in the L — > oo limit to be Gaussian- 
distributed around q^ v = 0, and so is c/diag- Now, let us 
hypothesize that there exists a spin-glass ordered state 
with a nonzero Edwards-Anderson SG order parameter 
9sg > 0- Reflecting the fact that (faiag transforms non- 
trivially under independent 0(2) rotations on the two 
replicas, even a self-overlap has nontrivial weights in 
^(<7diag) other than at ±ggQ. In the L — » oo limit, the 
self-overlap part of -P(<?diag) should be given by 



P(<?diag) = TT^diag) + 



2n (gf A^2 



JSC 



9di 



(5) 



The derivation is straightforward: When the two (essen- 
tially identical) ordered states in the two replicas are con- 
nected via a proper global rotation of the rotation-angle 
O, the diagonal overlap is given by ?diag = <?sg cos ©- 
Uniformly distributed then gives the second term of 
Eq.(5). When the two ordered states are connected via 
an improper global rotation, the diagonal overlap can be 

given by gdiag = 9sg J2i cos(2a — 28^)/N, where a is an 
angle of the reflection axis with respect to the x-axis in 
spin space, and 9^ denotes the direction of the i-th spin 
in replica 1. Since spins should be oriented randomly 
on long length scale in the SG ordered state, Qdiag given 
above should vanish in the L — > oo limit for arbitrary 
a, contributing a delta function at <7diag = 0, the first 
term of Eq.(5). If the SG orderd state hypotesized here 
accompanies RSB, the associated nontrivial contribution 
would be added to the one given by Eq.(5). In any case, 
an important observation here is that, as long as the or- 
dered state possesses a finite SG LRO, the diverging peak 
should arise in P((7diag) at (?diag = ±9sg > irrespective the 
occurrence of the RSB. 

We show in Fig. 3(b) the calculated P(gdiag) at J/T = 
2.95, well below X"cg- The calculated P(gdiag) exhibits 
symmetric shoulders at nonzero values of gdiag, but as 
shown in the inset, these shoulders get suppressed with 
increasing L, not showing a divergent behavior. We 
also calculate P(^diag) up to a still lower temperature 
J/T = 7, though only for smaller lattices L = 4, 6, 8, to 
observe an essentially similar behavior. Hence, we con- 
clude that the chiral-glass ordered state does not accom- 
pany the standard SG order with nonzero q§Q, at least 
up to temperatures » Xcg/3. No evidence of the suc- 
cessive chiral and spin transitions suggested in Ref. JTs|] 
was found. Strictly speaking, the observed suppression of 
the shoulder is still not inconsistent with the Kosterlitz- 
Thouless (KT)-like critical ordered state. However, such 
a critical SG ordered state is not supported by our data 
of g s in Fig. 1(b). Furthermore, the spin-glass exponent 
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reported in Ref. E^] assuming the simultaneous spin and 
chiral transition, v a = 1.2(4), is far from from the lower- 
critical-dimension value, v = oo, generically expected in 
such a KT transition. 
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FIG. 3. Chiral-overlap distribution function (a), and di- 
agonal spin-overlap distribution function (b), at J/T — 2.95 
well below Tqg- Inset of Fig.(b) is a magnified view of the 
shoulder part of P(gdiag). 

Finally, we wish to refer to the possible cause why si- 
multaneous spin and chiral orderings were apparently ob- 
served in certain numerical simulations. We believe this 
would closely be related to the length and time scale of 
measurements. We stress that the spin-chirality decou- 
pling is a long-scale phenomenon: At short scale, the chi- 
rality is never independent of the spin by its definition, 
roughly being its squared (k ~ S 2 ). Hence, the behavior 
of the spin-correlation related quantities, including the 
SG order parameter itself which is a summed correlation, 
might well reflect the critical singularity associated with 
the chirality up to certain length and time scale. If so, 
apparent, not true, "spin-glass exponents" at T = Tqq 
would be v' s ~ v K ~ 1.2 and r)' s ~ —0.4, the latter being 
derived from the short-scale relation, 1 + r] K ~ 2(1 + rj' s ). 
However, such a disguised criticality in the spin sector is 
only a short-scale phenomenon, not a true critical phe- 
nomenon. The length and time scales above which the 



spin-chirality separation becomes evident in correlations 
should roughly be given by the (finite) SG correlation 
length and correlation time at T = Tqg- We estimate 
these scales roughly of order 10 lattice spacings and some 
10 5 MC time steps (in the standard Metropolis dynam- 
ics). Meanwhile, the reason why the spin-chirality decou- 
pling looks evident already for smaller lattices in other 
types of quantities such as -P(<?diag) and g s remains to be 
understood. 

The numerical calculation was performed on the Hi- 
tachi SR8000 at the supercomputer center, ISSP, Univer- 
sity of Tokyo. MSL was supported by the Polish agency 
KBN (Grant no. 2P03B-146-18). 
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